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ABSTRACT 

Dark energy is an important science driver of many upcoming large-scale surveys. With small, 
stable seeing and low thermal infrared background, Dome A, Antarctica, offers a unique opportunity 
for shedding light on fundamental questions about the universe. We show that a deep, high-resolution 
imaging survey of 10,000 square degrees in ugrizyJH bands can provide competitive constraints on 
dark energy equation of state parameters using type la supernovae, baryon acoustic oscillations, and 
weak lensing techniques. Such a survey may be partially achieved with a coordinated effort of the 
Kunlun Dark Universe Survey Telescope (KDUST) in yJH bands over 5000-10,000 deg 2 and the 
Large Synoptic Survey Telescope in ugrizy bands over the same area. Moreover, the joint survey can 
take advantage of the high-resolution imaging at Dome A to further tighten the constraints on dark 
energy and to measure dark matter properties with strong lensing as well as galaxy-galaxy weak 
lensing. 

Subject headings: cosmological parameters — distance scale — gravitational lensing — large-scale 
structure of universe 
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1. INTRODUCTION 

Antarctic Plateau, especially the Kunlun Station un- 
der construction by the Polar Research Institute of China 
(PRIC), provides a unique opportunity for wide field as- 
tronomical surveys targeting cosmological studies. Astro- 
nomical site survey of Dome A, Antarctica was enabled 
by the International Polar Year (IPY) endorsed PANDA 
program (Yang ct al. 2009) led by the PRIC and the 
Chinese Center for Antarctica Astronomy (CCAA). Sev- 
eral international teams contributed to this effort. In par- 
ticular, the power and on-site laboratory system built by 
the University of New South Wales (BUN'S) has pro- 
vided the platform for all the site survey instruments. 

In its two years' operation, the site survey effort proves 
practically all aspects of the theoretical expectations of 
the Dome A site for astronomical observations. Prelim- 
inary analyses show that the boundary layer of atmo- 
spheric turbulence to be around 10 — 20 meters dur- 
ing the Antarctic winter (Ashley et al. 2010). Similar 
to the relatively better studied neighboring Dome C 
site (Fossat ct al. 2010), the Dome A site may enjoy free 
atmospheric seeing conditions of about 0.3 arcsec seeing 
above this boundary layer, thus making it an ideal site 
for high angular resolution wide area surveys. 

The temperature at Dome A is around —60 to — 70°C, 
making it the coldest spot on the surface of the Earth. 
This implies a very low thermal background emission in 
the thermal infrared. The Dome A site is thus also the 
best site for astronomical observations in the near in- 
frared wavelength region. 

One other especially exciting property of the site is the 



lack of water vapors due to the high altitude and the low 
temperature. This implies that the site is also ideal for 
terahertz observations which have been impossible from 
any temperate sites on Earth. 

While quantitative site properties are still under anal- 
yses and longer term monitoring are still needed to 
firmly establish the astronomical potential of the Dome 
A site, there is no doubt that a survey project at 
Dome A can be highly complimentary to programs 
such as the Large Synoptic Survey Telescope 1 (LSST, 
LSST Science Collaborations 2009), the Joint Dark En- 
ergy Mission 2 , and Euclid 3 . For example, a survey at 
Dome A can provide near infrared data that compliments 
the deep optical band survey of the LSST; a deep survey 
in the near infrared combined with the optical data from 
LSST can reveal high redshift objects at z ~ 10, which 
are not detectable in the optical. The Kunlun Dark Uni- 
verse Survey Telescope (KDUST) is a 6-to-8-meter wide- 
area survey telescope being designed by the CCAA. The 
preliminary design includes a 3 x 3 square degree opti- 
cal camera with 0".15 pixel, and an infrared camera of 
lxl square degree at 0".l /pixel optimized for 1 — 3.5 
/xm surveys. 

One of the key science missions of KDUST is to in- 
vestigate the mystery of the accelerated cosmic expan- 
sion (Riess et al. 1998; Perlmutter et al. 1999) using mul- 
tiple techniques. In this paper, we estimate how well 
an ideal 10,000 deg 2 ugrizyJH survey can constrain the 



1 Sec http://www.lsst .org/. 

2 See http://jdem.gsf cnasa.gov/. 

3 Sec http://sci.esa.int/euclid/. 
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dark energy equation of state (EOS) with weak lensing 
(WL), baryon acoustic oscillations (BAOs), and type la 
supernova (SNc) luminosity distances. These dark en- 
ergy probes have different sensitivities to the cosmic 
expansion and structure growth as well as various sys- 
tematic uncertainties in the observations, and hence are 
highly complementary to each other for constraining dark 
energy properties (e.g., Knox et al. 2005; Zhan 2006; 
Zhan et al. 2009). 

The SNe technique relies on the standardizable candle 
of the SNe intrinsic luminosity (Phillips 1993) to measure 
the luminosity distance, -Dl(z). Dark energy properties 
can then be inferred from the distance-redshift relation. 
The BAO technique utilizes the standard ruler of the 
baryon imprint on the matter (and hence galaxy) power 
spectrum (Peebles & Yu 1970; Bond & Efstathiou 1984) 
to measure the angular diameter distance, Dj^(z), and, if 
the redshifts are sufficiently accurate, the Hubble param- 
eter, H(z) (Eisenstein et al. 1998; Cooray et al. 2001; 
Blake & Glazebrook 2003; Hu & Haiman 2003; Linder 
2003; Seo & Eisenstein 2003). The WL technique has 
the advantage that it can measure both Da{z) from 
the lensing kernel and the growth factor of the large- 
scale structure G(z) (Hu & Tegmark 1999; Huterer 2002; 
Refregier 2003; Takada & Jain 2004; Knox et al. 2006a; 
Zhan et al. 2009). 

Because of the excellent seeing condition and infrared 
accessibility at Dome A, KDUST has a number of ad- 
vantages for the commonly used cosmological probes. For 
example, the signal-to-noisc ratio for point sources is in- 
versely proportional to the seeing. Thus, a 6 meter tele- 
scope at Dome A (~ 0".3 median seeing in the optical) 
would be equivalent to a 14 meter telescope at a temper- 
ate site (~ 0".7 seeing) for point-source observations with 
the same sky background level. An 8m KDUST could 
detect SNe out to redshift 3 in the -Kdark (2.27-2.45yum, 
redward of K) band (Kim et al. 2010). Although high-z 
distances are not sensitive to conventional dark energy, 
they can be used to determine the mean curvature accu- 
rately, which, in turn, helps constrain dark energy EOS 
(Linder 2005; Knox et al. 2006b). Moreover, even though 
dark energy is thought to be sub-dominant at high red- 
shift, there is no direct evidence to prove one way or an- 
other. Measurements of SNe at z > 2 will provide crucial 
data for tests of early dark energy. 

Small and stable seeing is particularly helpful for WL. 
One could resolve more galaxies at the same surface 
brightness limit, which reduces the shape noise for shear 
measurements. Fine resolution helps measure the shape 
accurately and reduce the shear measurement system- 
atic errors. In addition, deep JH photometry can track 
the 4000A break of an elliptical galaxy to z ~ 3 and 
improve photometric redshifts (photo-zs) as well as sys- 
tematic uncertainties in the photo-z error distribution 
(Abdalla et al. 2008), which has a large impact on WL 
constraints on the dark energy EOS (Huterer et al. 2006; 
Ma et al. 2006; Zhan 2006). 

Adding K or Kdark band will certainly improve galaxy 
photo-zs, especially at z > 3. However, currently planned 
multiband dark energy surveys use galaxies at z < 3, 
and their concern is the confusion between z < 0.5 el- 
lipticals and 2 < z < 3.5 star-forming galaxies, which 
is greatly mitigated by u and JH bands (Abdalla et al. 



2008). Another consideration is that Dome A is far more 
advantageous at -Kdark band than at K band because of 
the low thermal background there, so -Kdark is likely to 
be chosen over K. This would leave a considerable gap 
in the wavelength coverage and reduce the already-small 
gain on photo-zs in the useful redshift range for dark en- 
ergy investigations. Therefore, we do not discuss utilities 
of wavebands beyond H in this paper. Nevertheless, the 
-Kdark band is crucial for a broad range of other sciences 
and will be an important aspect of the KDUST survey. 
This paper is organized as follows. Section 2 dis- 
cusses the survey plan of KDUST including that of its 
pathfinder in context of the LSST survey. We then con- 
sider a joint KDUST and LSST survey in section 3 for 
constraining the dark energy EOS with BAO, WL, and 
type la SN techniques. The results are presented in both 
two-parameter space where the dark energy EOS is pa- 
rameterized as w(z) = wq + w a z(l + z) _1 and in model- 
independent principle component space. The conclusion 
is drawn in section 4. 

2. A POTENTIAL KDUST SURVEY PLAN 

Dome A has a great potential for cosmology, but, given 
other ambitious projects that will be concurrent with 
KDUST, one must carefully plan the KDUST survey to 
make the best use of the Dome A site. As we discuss, 
below, a potentially efficient strategy for KDUST is to 
focus on the near infrared (NIR) bands and incorporate 
optical data from LSST or other surveys. 

In the optical bands, a 6m KDUST is about twice as 
fast as LSST for surveying sky-dominated point sources 
at the same sky level, in which case the survey speed 
is proportional to the aperture and field of view and in- 
versely proportional to the seeing disk area. In reality, au- 
rorae increase the sky brightness in short wavelengths. It 
is estimated that the sky brightness at Dome A would be 
twice as bright as that at the best temperate site in B and 
20%-30% brighter in V (Saunders et al. 2009); measure- 
ments have shown that the median i-band sky brightness 
at Dome A in 2008 4 was 19.81 mag arcsec -2 and 20.46 
mag arcsec -2 during dark time, better than that at other 
good sites (Zou et al. 2010). With the above considera- 
tions, the 6m KDUST could survey 10,000 deg 2 to LSST 
depths in griz and much deeper in y (~ 26 mag, 5c point 
sources) in 2.5 years. The NIR camera of KDUST would 
have a much smaller field of view. It could reach J = 25 
mag and H = 24.6 mag over the same area in 3.5 years. 

Since LSST plans to survey the southern sky in ugrizy, 
it is not absolutely necessary for KDUST to survey in the 
optical again except in the y band. LSST would spend 
20% of its time in y band to achieve a 5-tr limiting mag- 
nitude of 24.4 for point sources, which is 2.8 magnitudes 
shallower than its r band limit. From photo- z consider- 
ation, it is desirable to have the y band limit not too 
much shallower than the limits in shorter wavebands. 
KDUST could improve the situation in its 10,000 deg 2 
survey area, which would be covered by LSST as well. 
A joint effort of KDUST in yJH and LSST in ugrizy 
would save both projects time while achieving better per- 
formance in measuring the dark energy EOS. 

Here, we use SNAP and LSST as references to model 

4 The sky background was seriously affected by the Moon in 
2008, as it was always close to full when above the horizon. 
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Table 1 

Model parameters of survey data and priors. 





Ideal 10k 


LSST 


KDUST+LSST 3 


area (deg 2 ) 


10,000 


20,000 


10,000 + 10,000 


rig (arcmin -2 ) 


70 


40 


60 40 


2* 


0.5 


0.3 


0.4 0.3 


TzO 


0.03 


0.05 


0.03 0.04 


a P (6z) 


0.2ct 2 


0.3cr z 


0.2cr z 0.2cr z 


M/7) 


0.002 


0.005 


0.002 0.003 


(A7) 2 (xlO- 10 ) 


4 


10 


4 6 



a We assume that KDUST would survey 10,000 deg 2 in yJH bands 
and several hundred deg 2 in griz bands around the south celestial 
pole, excluding areas with heavy galactic extinction. In the 10,000 
deg 2 overlap region between KDUST and LSST, the combined sur- 
vey is close to the ideal 10k survey. The data outside the overlap 
region come solely from LSST but are expected to have smaller 
systematics than they do without KDUST because of better cali- 
brations from the overlap region. 

an ideal 10,000 square degree survey (named "ideal 10k") 
in ugrizyJH bands to comparable depth as SNAP. We 
assume that its galaxy redshift distribution follows 



n(z) ex z e 



z/z 



with 



= 0.5. The projected galaxy number density 
n g is 70 arcmin -2 , and the distribution peaks at 2z*. 
For LSST, we adopt z* = 0.3 and n g = 40 arcmin -2 
(LSST Science Collaborations 2009). The model param- 
eters of the survey data are summarized in Table 1. 

The galaxy distribution rii(z) in the ith bin is sampled 
from n(z) by (Ma et al. 2006; Zhan 2006) 



rii(z) = n(z)V{z pi ,z pl ; 



•-)• 



where the subscript p denotes photo- z space, z pi and z pi 
define the extent of bin i, and Via, b; z) is the probability 
of assigning a galaxy that is at true redshift z to the 
photo-z bin between z v = a and b. We approximate the 
photo-z error to be Gaussian with bias 5z and rms a z = 
020 (1 + •?), and the probability becomes 



'P^p.i, Zp.il z) = I(Zp,i, Zp/, z)/I(0, 00; z), 



I(a,b;z)-- 




dz p exp 



(z p — z — 5zY 



2a? 



The normalization 1(0, 00; z) implies that galaxies with a 
negative photo-z have been excluded from n(z). We use 
40 photo-z bias Sz and 40 photo-z rms a z parameters 
evenly spaced between z = and 5 to model the photo- 
z error distribution in the z-z p space; the photo-z bias 
and rms at any redshift are linearly interpolated from 
these 80 parameters. Note that the photo-z parameters 
are assigned independent of galaxy bins. We assume that 
a z0 = 0.03 for the ideal survey and a z o = 0.05 for LSST. 
Uncertainties of the photo-z parameters (or, the photo- 
z error distribution in general) have a large impact on 
the dark energy constraints from WL and are referred 
to as photo-z systematics. Therefore, the prior on the 
photo-z error distribution is an important quantity to 
specify when reporting WL constraints on dark energy. 
To reduce the dimension of the investigation, we set a 
simple function for all the priors on the photo-z bias 
parameters, ap(dz) — 0.2a z , and peg the priors on 
photo-z rms parameters to those on the bias parameters: 



crp(cr z ) = v2ap(5z). For Gaussian photo-z errors, these 
priors correspond to a calibration sample of 25 spectra 
per redshift interval of the photo-z parameters. We set 
ap(<5z) = 0.3er z , reflecting larger uncertainties with fewer 
filter bands. However, such difference in the photo-z pri- 
ors is not important if one performs a joint analysis of 
the BAO and WL techniques, which takes advantage 
of the self-calibration of photo-z error distribution by 
galaxy power spectra (Schneider et al. 2006; Zhan 2006; 
Zhang et al. 2009, referenced herein) 

LSST will benefit from KDUST yJH data in a num- 
ber of ways. In the 10,000 deg 2 overlap region, KDUST 
data will (1) improve photo-zs directly, (2) increase the 
galaxy sample with those that would not meet the LSST 
optical photometry selection criteria, or whose photo-z 
could not be reliably determined in the absence of the 
deep yJH data, or whose shape could not be measured 
well because of the larger seeing at the LSST site, and 
(3) provide calibration of various systematic errors such 
as those in shear measurements. Even in the non-overlap 
region, LSST could still improve photo-zs and shear mea- 
surements because of the calibration within the overlap 
region. This is reflected in the last column of Table 1 
where we construct a joint survey combining LSST and 
KDUST. 

For the SN data, we assume that KDUST would ob- 
tain a SNAP-like sample of 2000 SNc reaching redshift 
1.7 as well as a sample of 1000 local and nearby SNe. De- 
tails of the assumptions for the SNAP experiments and 
the calculation of the Fisher matrices can be found in 
(Pogosian et al. 2005) and in (Zhao et al. 2009). 

3. CONSTRAINING DARK ENERGY WITH KDUST AND 

LSST 

In this section, we estimate the constraints on the dark 
energy EOS from the joint survey of KDUST and LSST 
using BAO, WL, and SN techniques. 

3.1. Baryon Acoustic Oscillations and Weak Lensing 

The angular power spectra of the galaxy number den- 
sity n(9) and (E-mode) shear j(0) can be written as 
(Hu & Jain 2004; Zhan 2006) 

P # K W = 5 / dzH(z)D A (z)W?(z)Wj(z)tf s (k;z), 

(1) 
where lower case subscripts correspond to the tomo- 
graphic bins, upper case superscripts label the observ- 
ables, e.g., X = g for galaxies or 7 for shear, A 2 (/c;z) 
is the dimensionlcss power spectrum of the density field, 
and k = £/D A (z). The window functions are 



Wf(z)=b(z) 



m(z) 



W?(z) = - 



3 n m # 2 D A (z) 



2 H(z) 



a c 



dz 



, nj( z') D A {z,z') 
H D A {z>) ' 



where b(z) is the linear galaxy clustering bias, and Q m 
and Hq are, respectively, the matter fraction at z = and 
Hubble constant. The galaxy redshift distribution n^z) 
in the ith tomographic bin is an average of the under- 
lying three-dimensional galaxy distribution over angles, 
and the mean surface density fii is the total number of 
galaxies per stcradian in bin i. For WL, we use 10 bins 
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Figure 1. Left panel: Five galaxy auto power spectra from the ideal 10k survey. The BAO feature is visible around £ ~ 400 in the 
highest redshift bin and move toward lower multipoles in lower redshift bins. Note, however, that the linear scale of the BAOs is fixed in 
the comoving frame. The amplitude of the power spectra and that of the BAO feature decreases as the photo- z rms error increases. The 
shaded area gives the 1-<t error (e.g., sample variance and shot noise) of the power spectrum per £ in the bin centered at z = 1.66. The 
power spectra are shifted for clarity. Right panel: Cross power spectra between the bin at z = 1.66 and its second neighboring bins and 4th 
neighboring bins. 
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Figure 2. Same as Figure 1 but for the shear power spectra in 5 different redshift bins. The shear power spectra are smoother than the 
galaxy power spectra, because lensing kernels are much broader that galaxy distributions in the photo-z bins. 



evenly spaced between z = to 3.5, and for BAO, we use 
30 bins from z = 0.15 to 3.5 with bin width proportional 
to 1 + z. 

The lcr statistical error of the mean power spectrum 
in the multipole range (£, t + A£) is given by 



AP$ X (£) 



pirw> 



- z p xY (f) 

'Y (2^+l)A£/ sky « l h 



where 5^ h is the Kronecker delta function, / s k y is the sky 
coverage, g rms = 1, and 7 rms = 0.18 + 0.042z. 

Figure 1 shows several examples of the galaxy auto and 
cross power spectra, and Figure 2 shows shear auto and 
cross power spectra. Cross power spectra between fore- 
ground galaxy bins and background shear bins, i"S , are 
not shown, but they have similar characteristic shapes 
of those in Figure 1 and Figure 2. The amplitude of the 
galaxy power spectra and the strength of the BAO fea- 
ture decreases with increasing photo-z rms error, because 
a larger photo-z rms means more smoothing in the ra- 
dial direction. The amplitude of the shear power spec- 



trum increases as redshift increases, because higher red- 
shift galaxies are lensed by more intervening matter and, 
hence, have stronger shear signal fluctuations. The am- 
plitude of the galaxy cross power spectrum between two 
bins is very sensitive to the separation between the two 
bins in true-redshift space. Hence, galaxy cross power 
spectra can be used to calibrate the photo-z error distri- 
bution. 

3.2. Constraints on wq and w a 

We use the Fisher information matrix (Tegmark 1997) 
to estimate the errors of the parameters of interest. In 
summary, the Fisher matrix of the parameter set {q a } is 
given by 



F a/ 3 = / sk y 2_, 



2£+l „-!dC e ^-xdCi 



-TrC 



dq a 



dqp 



(2) 



AT 
ij 



with (Ci) 

imum marginalized error of q a is o~(q a ) = (F~ l )J a . In- 
dependent Fisher matrices are additive; a prior on q a , 



PSj (£) for galaxies and shear. The min- 

1/2 



ap(q a ), can be introduced via F^ 



F n 



\q a ) 



We extend the additive and multiplicative shear power 
spectrum errors in Huterer et al. (2006) to include the 
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galaxy power spectrum errors: 

{Cf Y ). l3 = (1 + 5%f x + 5%fJ)P^{e) + 



>XY 



cK^rms X aX aY 



,(3) 



where p x determines how strongly the additive errors of 
two different bins are correlated, and -q x and £ x account 
for the scale dependence of the additive errors. Note that 
the multiplicative error of galaxy number density is de- 
generate with the galaxy clustering bias and is hence ab- 
sorbed by hi . At the levels of systematics future surveys 
aim to achieve, the most important aspect of the (shear) 
additive error is its amplitude (Hutcrcr ct al. 2006), so 
we simply fix p x = 1 and r/ x = 0. For more compre- 
hensive accounts of the above systematic uncertainties, 
see Huterer et al. (2006); Jain et al. (2006); Ma et al. 
(2006); Zhan (2006). 

Forecasts on dark energy constraints are sensitive to 
the priors on the shear multiplicative errors (cp(/7)) and 
the amplitudes of the shear additive errors (A]), so we 
list them in Table 1. See Wittman (2005); Massey ct al. 
(2007); Paulin-Hcnriksson et al. (2008) for detailed work 
on these systematic uncertainties and Zhan et al. (2009) 
for a discussion on erp(/ 4 7 ) and A] for LSST. It has 
been demonstrated that small and stable point spread 
function in space leads to better shear measurements 
(Kasliwal et al. 2008). In the absence of a detailed in- 
vestigation, we simply choose values or priors of these 
systematic error parameters somewhat arbitrarily be- 
tween what might be achieved by space projects and 
what have been used in LSST forecasts. We infer from the 
Sloan Digital Sky Survey galaxy angular power spectrum 
(Tcgmark et al. 2002) that the additive galaxy power 
spectrum error due to extinction, photometry calibra- 
tion, and seeing will be at (Af) 2 = 10~ 8 level. Since the 
amplitude of Af is fairly low compared to the galaxy 
power spectra, its value does not affect the forecasts 
much. 

In summary, the parameter set includes 11 cosmologi- 
cal parameters and 170 nuisance parameters. The cosmo- 
logical parameters are wq, w a , the matter density ui m , the 
baryon density Wb , the angular size of the sound horizon 
at the last scattering surface 9 S , the curvature parame- 
ter £7k, the scalar spectral index n s , the running of the 
spectral index a s , the primordial Helium fraction Y p , the 
election optical depth r, and the normalization of the 
primordial curvature power spectrum Ap. Note that Y p 
and r are solely constrained by the cosmic microwave 
background (CMB), which is introduced as priors. The 
nuisance parameters include 40 photo- z bias parameters, 
40 photo- z rms parameters, 40 galaxy clustering bias pa- 
rameters, 30 galaxy additive noise parameters, 10 shear 
additive noise parameters, and 10 parameters for shear 
calibration errors. We use multipoles 40 < £ < 2000 for 
WL and 40 < I < 3000 for BAO. In addition, we require 
A 2 (£/Da;z) < 0.4 for BAO to reduce the influence of 
nonlinear evolution. The lower cut in £ is to minimize 
the dependence of the forecasts on particular models of 
dark energy perturbation and the integrated Sachs- Wolfe 
effect, which affect only very large scales. 

Figure 3 presents the forecasts of la error contours 
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Figure 3. Forecasts of lex errors on the dark energy EOS pa- 
rameters wq and Wo, for the ideal 10k survey WL (solid line), 
BAOs (dotted line), SNe (dashed line), and the three combined 
(shaded area). We have included Planck priors in all the results. 
Although the CMB priors have a significant impact on the SN re- 
sults and to a lesser degree on WL results and BAO results, they 
have a much smaller effect on the joint constraints of BAO+WL 
and BAO+WL+SNe. 
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Figure 4. Forecasts of 1<t errors on the dark energy EOS param- 
eters wo and U) a for LSST using BAO+WL (dash-dotted line), the 
combination of LSST and KDUST using BAO+WL (solid line, see 
Table 1), the combination of LSST and half of KDUST (labeled as 
KDUST5k) using BAO+WL (dotted line), and the combination of 
LSST and KDUST5k using BAO+WL+SNe (shaded area). 

of the dark energy EOS parameters w$ and w a for the 
ideal 10k survey WL (solid line), BAO (dotted line), SNe 
(dot-dashed line), joint BAO and WL (dashed line), and 
the three combined (shaded area). Since different probes 
have different parameter degeneracy directions in the full 
parameter space, including both cosmological parameters 
and nuisance parameters (such as the photo-z bias and 
rms parameters), a joint analysis can reduce the error 
significantly. One example is that the BAO technique 
can determine the curvature parameter and the matter 
density far better than the SNe technique, whereas the 
latter needs strong priors on fix and w m in order to place 
tight constraints on wq and w a (Lindcr 2005; Knox ct al. 
2006b). Another example is that the WL technique is 
sensitive to the priors on the photo-z parameters, which 
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can be calibrated by the BAO technique. By comparing 
the results in Figure 3 with the LSST BAO+WL result 
in Figure 4, one sees that the ideal 10,000 deg 2 ugrizyJH 
survey can place comparable constraints on the dark en- 
ergy EOS parameters to the LSST survey. 

We show the dark energy constraints from combina- 
tions of LSST with the KDUST 10,000 deg 2 survey 
and LSST with half of the KDUST survey (labeled as 
KDUST5k) in Figure 4. The KDUST5k survey could be 
carried out by a 2.5m KDUST pathfinder in 12 years. We 
assume that 5000 deg 2 high-resolution yJH imaging is 
enough to reach the systematic calibration floor, so the 
improvement to the LSST survey outside the KDUST 
survey area is kept the same for both proposals. Un- 
der this condition, the WL+BAO constraints on wq and 
w a from LSST+KDUST5k are nearly the same as those 
from LSST+KDUST; they improve the LSST-alone dark 
energy task force (DETF, Albrecht et al. 2006) figure 
of merit (FOM) by 30%. Adding SNAP-like SN data 
from KDUST as well can significantly increase the DETF 
FOM; the result is as good as LSST or LSST+KDUST 
joint BAO and WL constraints on wq and w a in the ab- 
sence of systematics. 

The current constraint on {wo, w a } is (Zhao & Zhang 
2010), 



Wq 



n Qn +0. 11+0.23 
u ' JU -0.11-0.22' 



n 9/1+0.56+0.98 
" u - z ^-0.55-1.2 • 



(4) 



derived from a joint Markov Chain Monte Carlo 
(MCMC) analysis of the recently released SNe "Con- 
stitution" sample (Hicken et al 2009), the WMAP five- 
year data (Komatsu et al 2009) 5 , and the Sloan 
Digital Sky Survey (SDSS) Luminous Red Galaxy 
(LRG) sample (Tegmark et al 2006). In this analy- 
sis, the dark energy perturbation (DEP), which is 
important in the parameter estimation (Zhao et al. 
2005; Fang, Hu & Lewis 2008), was consistently in- 
cluded based on the treatment developed in (Zhao et al. 
2005). The central values indicate that the 'quintom' 
scenario (Feng, Wang & Zhang 2005) is mildly favored, 
namely, the EOS today w(z)\ z=0 = w > — 1, while EOS 
in the far past w(z)\ z=0 o = wq + w a < — 1. This is consis- 
tent with the recent published result using the 'Constitu- 
tion' SNe sample (Shafieloo, Sahni & Starobinsky 2009; 
Biswas & Wandelt 2009; Qi, Lu & Wang 2009; Wei 
2009; Huang et al. 2009). The error bars in Equation 4 
can be used to estimate the FOM of current data, and we 
find that KDUST+-LSST can improve the current FOM 
by two orders of magnitude. 

3.3. Principal Component Analysis on w(z) 

To investigate the constraints on dark energy 
equation-of-state w{z) from the surveys in a model- 
independent way, we follow (Huterer & Starkman 2002; 
Crittenden, Pogosian & Zhao 2009) and employ a Prin- 
ciple Component Analysis (PCA) approach. 

We choose 40 uniform redshift bins, stretching to a 
maximum redshift of z — 3, and allow the high redshift 
(z > 3) equation of state to vary. To avoid w(z) with in- 
finite derivatives, each bin rises and falls following a hy- 
perbolic tangent function with a typical transition width 

5 We have checked that the result is largely unchanged if the 
WMAP 5 year data is replaced with WMAP 7 year data. 



dz of order 10% of the width of a bin. We choose a con- 
stant w = — 1.0 as the fiducial model, which is consistent 
with all the present data. Since dark energy perturba- 
tions play a crucial role in the parameter estimation, we 
use a modified version of CAMB which allows us to calcu- 
late DEP for an arbitrary w(z) consistently (Zhao et al. 
2005). 

For the PCA, we calculate the Fisher matrices based on 
four kinds of observables: supernovae, CMB anisotropics, 
galaxy number counts (GC) correlation functions and 
WL observations. We also include all the possible cross- 
correlations among these, including CMB x galaxy and 
CMBxWL, which are sensitive to the integrated Sachs- 
Wolfe effect, as well as galaxy- weak lensing correlations. 

We first calculate the Fisher matrices for each of the 
observables, F^, where the indices i,j run over the pa- 
rameters of the theory, in our case the binned Wi(z). We 
then find the normalized eigenvectors and eigenvalues of 
this matrix {ei(z), Aj}, and write 



F = W T AW, 



(5) 



where the rows of W are the eigenvectors and A is a 
diagonal matrix with elements A^. The Fisher matrix is 
an estimate of the inverse covariance matrix we expect 
the data to give us and the eigenvalues reflect how well 
the amplitude of each eigenvector can be measured. The 
true behavior of the equation of state can be expanded 
in the eigenvectors as 



/v 



w(z) = U> fid (z) + ^ a * e i( z ) 



(6) 
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and the expected error in the recovered amplitudes is 

given by <x(a,) = A ; . 

For our forecasts, we assume three surveys: the ideal 
10k survey, LSST and KDUST+LSST as listed in 
Table 1. We also combine the Planck survey for CMB and 
SNAP SNe. We assume a flat universe and marginalize 
over the intrinsic SN magnitude M and the galaxy bias 
parameters. 

In the upper panel of Figure 5, we show the spectra 
of eigenvalues a~ 2 (ai) of the Fisher matrices for these 
three tomographic surveys, and in the lower panel, we 
show the ratio of the eigenvalues for LSST+KDUST to 
that for LSST. As we can see, the ideal 10k survey is 
competitive to LSST on the dark energy constraints, and 
adding KDUST to LSST can improve the constraints on 
the first few eigenmodes by as much as 85%. The first 
five best determined eigenvectors are shown in Figure 6. 
As we can see, the Nth eigenmode has N — 1 nodes in z, 
and the 'sweet spot' - the redshift where the uncertainty 
of w(z) gets minimized - is at z ~ 0.2. This is consistent 
with the analysis done in (Huterer & Starkman 2002; 
Crittenden, Pogosian & Zhao 2009). We also find that 
adding KDUST to LSST makes the eigenmodes stretch 
to slightly higher z because of the higher redshift reach 
with the addition of KDUST data. 

As stated in (Crittenden, Pogosian & Zhao 2009), all 
of the eigenvectors are informative, no matter how large 
the error bars are, if we have no prior knowledge of the 
possible w(z) behaviors. However, even without a phys- 
ical model for w(z), we would still be surprised if w{z) 
were much too positive (w ^> 1/3) or much too nega- 
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Figure 5. Upper panel: The eigenvalues (a~ 2 (ai)) for the raw 
Fisher matrices (no priors assumed) for ideal 10k survey, LSST, and 
LSST+KDUST. The grey shaded region shows the diagonal prior 
of cr m > 0.3; Lower panel: The ratio of the eigenvalues (a~ (q;)) 
for LSST+KDUST to that for LSST. 




Figure 6. The first five best determined eigenvectors for the ideal 
10k survey (dashed lines), LSST (solid lines), and LSST+KDUST 
(dash-dotted lines). No priors have been assumed and the ampli- 
tudes are normalized to unity. 

tive (w <C —1). This motivates us to choose some the- 
oretical priors to roughly separate the eigenmodes into 
those which are informative relative to the priors and 
those that are not. To apply the theoretical prior, we 
follow (Crittenden, Pogosian & Zhao 2009) to choose a 
correlation function describing fluctuations of w(z) away 
from the fiducial model. This smoothness prior can filter 
out the high frequency modes, while the low frequency 
modes remain unaffected. Also, as long as there are suf- 
ficient bins compared to the correlation length, the prior 
largely wipes out dependence on the precise choice of 
binning. 
As elaborated in (Crittenden, Pogosian & Zhao 



2009), the deviations of the equation of state from its 
fiducial model can be encapsulated in a correlation 
function: 

£ w (\z - z'\) = ((w(z) - w Rd (z))(w(z') - itffidCO)) • (7) 

and the covariance matrix of the binned equation of state 
is then 



(Si 



i §W i) = ^2 



Zi+A 



dz 



Zj+A 



dz'^ w {\z~z'\). (8) 



where the i bin is from z% to Zi + A, and we assume 
that all bins have the same width A = Zi+i — Zi. The 
variance of the mean equation of state over all the bins 
is, 

■" dz' 



dz 



-£w{z 



(9) 



And we assume that £ w (z) = £„,(0)/(l + (z/z c ) 2 ). where 
£io(0) is the variance of w at any given point and z c is 
the correlation length. 

To apply the prior, we need to specify the correlation 
length z c and tune £«i(0) so that the error in the mean, 
a m is in the range of [0.2,0.5] in order to be consistent 
with the observational uncertainty. We tried correlation 
lengths in the range 0.1 < z c < 0.4, where the upper 
limit was beginning to be strong enough to impact the 
observed modes for the SN. Then the resulting prior takes 
the form, 



oc exp 



1 



H 



wf d )CrM w 



10 - wf ) 



(10) 



where CV,- = (SwiSwj). This prior naturally constrains 
the high frequency modes without over constraining the 
lower frequency modes that are typically probed by ex- 
periments. 

Note that assuming no correlations among the bins 
is equivalent to using a delta function for the corre- 
lation prior, e.g. £(z) = £o^(z)- In such a case, one 
finds (SwiSwj) — £o<%/A and the mean variance is 
°m = £o/ z max- Thus, assuming a fixed total range, 
the bin variance should grow with the number of bins 



(ft 



CjjjiVbms to keep the mean variance unchanged. 



In order to evaluate the value added by a given survey, 
we use the mean squared error (MSE) as a FOM. In the 
absence of any priors, this is expected to be 



MSE = Tr F~ 



(11) 



Taking into account priors, F in this expression is re- 
placed by C" 1 +F. The MSE arising from the prior alone 
is MSE = Tr C, and it is independent of the shape of the 
correlation function. Adding more experimental data re- 
duces the MSE, and the amount by which the MSE is 
reduced can be seen as a measure of how informative the 
experiment is. 

In Table 2, we vary £ w (0) and z c while holding the prior 
constraint on the mean fixed, a nl = 0.3. For the case of 
a diagonal prior, the MSE for the prior alone is given 
by N hlns x (Sw 2 ) = A^ ins cr^ = 144. Adding the data, 
one could significantly improve the constraints on some 
of the eigenmodes, thus reduce the MSE. As the prior is 
diagonal, the new eigenmodes are also eigenmodes of the 
prior with the same eigenvalue, given by Nbins&m = 3.6 
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Table 2 

The mean squared error for various priors and experiments. 



diagonal 



0.1 z c =0A 



no data 



144.0 



35.0 



11.5 



Ideal 10k 

LSST 

KDUST+LSST 



116.0 [7] 
112.0 [9] 
108.3 [10] 



21.2 
19.9 
18.1 



1.7 
4.5 

3.7 



Note. — The mean squared error is related to the number of 
well constrained modes. The priors arc normalized so that cr m = 0.3 
and the "diagonal" prior means no correlations exist between bins. 
For this diagonal prior, we give in brackets the inferred number of 
modes meaningfully constrained by the observations. 

(shaded region in the upper panels of Figure 5). The 
reduction of the MSE thus roughly tells us how many 
modes can be constrained compared to the prior. For ex- 
ample, adding LSST reduces the MSE from 144 to 112, 
meaning that there are (144 — 112)/3.6 ~ 9 modes can 
be constrained by LSST, and similarly, the ideal 10k sur- 
vey is able to constrain 7 modes. With KDUST+LSST, 
one can actually constrain 10 eigenmodes. These num- 
bers can also be counted in the upper panel of Figure 5 
above the prior threshold shown in shade. 

For a correlated prior, we can analogously put a lower 
bound on the number of modes constrained by data us- 
ing the MSE shown in Table 2. For instance, for the 
case of z c = 0.1 and z c = 0.4, LSST combined with 
KDUST can at least constrain (35 — 18.1)/3.6 ~ 4 and 
(11.5 — 3.7)/3.6 ~ 2 modes, respectively. As we see, the 
number of new modes estimated in this way is reduced 
as the prior correlation length is increased. This is as ex- 
pected - as z c increases, more high frequency modes will 
be filtered out by the smoothness prior. 

4. CONCLUSION 

Dome A offers a very competitive site for studying dark 
energy Given the amount of resources required to build 
a large telescope and run a massive survey there, one 
must give the highest priority to programs that cannot 
be easily carried out elsewhere. Thus, a reasonable strat- 
egy is to focus on NIR imaging and collaborate with other 
surveys for optical data. Using LSST as an example, we 
show that a high-resolution 5000-10,000 deg 2 KDUST 
survey in yJH bands could improve LSST BAO+WL 
constraints on the dark energy EOS parameters wq and 
w a by reducing the photo- z and shear measurement sys- 
tematics. A SNAP-likc SN sample plus a large local and 
nearby SN sample from KDUST would further boost the 
DETF FOM by more than a factor of two. 

In addition to forecasts for the WQ-w a paramctriza- 
tion, we also apply a PCA approach to investigate the 
constraints on the dark energy EOS w(z) in a model- 
independent way. We find that regarding the number of 
the constrained eigenmodes oiw(z), an ideal 10,000 deg 2 
ugrizyJH survey, combined with Planck, can constrain 7 
eigenmodes, while KDUST+LSST can allow us to con- 
strain 3 more modes. 

We have not discussed dark energy probes such as 
strong lensing, cluster counting, and higher-order statis- 
tics of the same galaxy and shear data, which could fur- 
ther tighten the constraints on the dark energy EOS. 
Strong lensing constrains dark energy through the time 
delay effect as well as counting of strong lenses. It is 
also an excellent probe of dark matter halo structures 



and, hence, can be used to measure dark matter particle 
properties. With high-resolution imaging, one could ex- 
tract more cosmological information from strong lensing 
observations. Therefore, Dome A could be particularly 
advantageous for strong lensing studies. 

Dark energy forecasts depend crucially on the assumed 
properties of the survey data, including all the systemat- 
ics. Dome A has many advantages over other ground sites 
and has an environment close to that in space. Hence, we 
use well-studied LSST and SNAP as references to make 
crude estimates of the data for this investigation. Further 
work and detailed modeling are needed to give a more re- 
alistic assessment of the Dome A site for studying dark 
energy 
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